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Abstract 



0^ ', Extended version of lectures given at the University of Karlsruhe and at Calc-2000 

' school in Dubna. Properties of HQET as a field theory, methods of calculation of 

HQET diagrams and some simple applications are explained in detail. These lectures 
, can be used as an additional chapter with any modern QCD textbook. 

o 

g : 1 HQET Lagrangian 

§ ■ Let's consider QCD with a heavy flavour Q with mass m and a number of light flavours. 
"^H , We shall be interested in problems with a single heavy quark staying approximately at 
Q-i! rest. More exactly, let uj <^ m be the characteristic momentum scale. We shall assume 
Q^l that the heavy quark has the momentum \p\ ^ uj and the energy \pq — m\ ^ uj; light 
quarks and gluons have momenta \ki\ < and energies \koi\ ^ uj. Heavy quark effective 
theory (HQET) is an effective fleld theory constructed to reproduce QCD results for such 
problems expanded up to some order in uj/m. In practice, only a few flrst orders in 1/m 
expansion are considered, because complexity of the theory grows fast with the order. 



5^ I Let's start from the QCD Lagrangian 

L = Q{tip-m)Q + --- (LI) 

where Q is the heavy quark fleld, and dots mean all the terms with light quarks and gluons. 
The free heavy quark Lagrangian Q{ip — m)Q gives the dependence of the energy on the 
momentum po = + p^. If we assume that characteristic momenta |p | -Cm, then we 
can simplify the dispersion law to po = It corresponds to the Lagrangian Q{i'~^QdQ—m)Q. 
In our class of problems, the lowest-energy state ("vacuum") consists of a single particle 
— the heavy quark at rest. Therefore, it is convenient to use the energy of this state m as 
a new zero level. This means that instead of the true energy po of the heavy quark (or any 
state containing this quark) we shall use the residual energy po = Po — m. Then the on-shell 
heavy quark has the energy po = independently on the momentum. The free Lagrangian 
giving such a dispersion law is QijodoQ- The spin of the heavy quark at rest can be 
described by a 2-component spinor Q (we can also consider it as a 4-component spinor 
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with the vanishing lower components: 70Q = Q). Reintroducing the interaction with the 
gluon field by requirement of the gauge invariance, we arrive at the HQET Lagrangian 

L = Q+iDoQ + --- (1.2) 

where all light-field parts (denoted by dots) are exactly the same as in QCD. The field 
theory ( |1.2| ) is not Lorentz-invariant, because the heavy quark defines a selected frame — 
its rest frame. 

The Lagrangian (p.. 21) gives the static quark propagator 



S{P) = ^^—^, S{x) = S{xo)S{x), S{t) = -te{t). (1.3) 

In the momentum space it depends only on po but not on p, because we have neglected 
the kinetic energy. Therefore, in the coordinate space the static quark does not move. The 



unit 2x2 matrix is assumed in the propagator ( |1.3|) . It is often convenient to use it as 



a 4 X 4 matrix; in such a case the projector onto the upper components is assumed. 
The static quark interacts only with Aq; the vertex is ig6f^t°'. 

The static quark propagator in a gluon field is given by the straight Wilson line 

S{x) = -i9{xo)6{x)Pexpig J A^dx^ . (1.4) 

Many properties of HQET were first derived in the course of investigation of renormaliza- 
tion of Wilson lines in QCD. In fact, the HQET Lagrangian was used as a technical device 
for investigation of Wilson lines. 

Loops of a static quark vanish, because it propagates only forward in time. In other 
words, in the momentum space, all poles in the po plane lie in the lower half-plane; closing 
integration contours upwards, we get zero. 

The Lagrangian ( p..2[) can be rewritten in covariant notations P, 

L„ = Q> ■ DQ^ + ■■■ (1.5) 

where the static quark field is a 4-component spinor obeying the relation = Qv, 
and is the quark velocity. The momentum p of the heavy quark (or any state containing 
it) is related to the residual momentum p by 

p = mv+p, |p'"|<^m. (1-6) 

The static quark propagator is 

1 + ^ 1 

S{P}= /-^ (1.7) 
2 p ■ V + t(} 

and the vertex is igvH"'. 



2 



One can watch [0 how expressions for QCD diagrams tend to the corresponding HQET 
expressions in the hmit m oo. The QCD heavy quark propagator is 

S(p) ^ = '"'1+ ^^[t! = ^IJ- + o(l). (1.8) 

pz _ zmp-v+p'^ 2 p-v \m J 

A vertex ig'jH"- sandwiched between two projectors may be replaced by igvH"- (one 
may insert the projectors at external heavy quark legs, too). Therefore, any tree QCD 
diagram equals the corresponding HQET one up to 0{p/m) terms. In loops, momenta can 
be arbitrarily large, and the relation ( |1.8|) can break. But if we renormalize HQET, loop 
integrals become convergent. In them, characteristic loop momenta are of order u, and 
one may use ( p..8|) . 

Renormalization properties (anomalous dimensions etc.) of HQET differ from those of 
QCD. The ultraviolet behavior of an HQET diagram is determined by the region of loop 
momenta much larger than the characteristic momentum scale of the process u, but much 
less than the heavy quark mass m (which tends to infinity from the very beginning). It has 
nothing to do with the ultraviolet behavior of the corresponding QCD diagram with the 
heavy quark line, which is determined by the region of loop momenta much larger than m. 
In the conventional QCD, the first region produces hybrid logarithms which are difficult to 
sum up. In HQET, hybrid logarithms become ultraviolet logarithmic divergences governed 
by the renormalization group with corresponding anomalous dimensions. 

The HQET Lagrangian ( |1.2| ) possesses the SU(2) heavy quark spin symmetry The 
heavy quark spin does not interact with gluon field at m — > oo, because its chromomagnetic 
moment is of the order of 1/m by dimensionality. Therefore, we can rotate the heavy quark 
spin separately. 

Not only the orientation, but also the magnitude of the heavy quark spin is irrelevant 
in HQET. We can switch off the heavy quark spin, making it spinless. This often simplifies 
counting form factors in HQET. Or, if we want, we can make the heavy quark to have spin 
1. This leads to a supersymmetry group called the superflavour symmetry It allows 
one to predict properties of hadrons with a scalar or vector heavy quark (appearing in 
some extensions of the Standard Model). This idea can also be applied to baryons with 
two heavy quarks, because they form a small size (of order l/{mas)) bound state with spin 
or 1, antitriplet in colour. 

HQET has great advantages over QCD in lattice simulation of heavy quark problems. 
Indeed, the applicability conditions of the lattice approximation to problems with light 
hadrons require that the lattice spacing is much less than the characteristic hadron size, 
and the total lattice length is much larger than this size. For simulation of QCD with a 
heavy quark, the lattice spacing must be much less than the heavy quark Compton wave- 
length 1/m. For b quark, this is technically impossible at present. The HQET Lagrangian 
does not involve the heavy quark mass m, and the applicability conditions of the lattice 
approximation to HQET are the same in the case of light hadrons. 
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2 One-loop massless propagator diagrams 



The main aim of these lectures is to discuss the methods used for calculation of Feynman 
diagrams in HQET, and some of the simplest applications. These methods have much in 
common with the ones used to calculate multiloop diagrams in massless theories. There- 
fore, for convenience of the readers, we shall recall some well-known facts about massless 
diagrams, too. 

In this Section, we shall calculate the one-loop massless propagator diagram with arbi- 
trary degrees of the denominators (Fig. [l|) 

where c? = 4 — 2e is the space-time dimension. 

k 




k + p 



Figure 1: One- loop massless propagator diagram 

The first step is to combine the denominators together. To this end, we write l/a" for 
Rea > as 



oo 

1 1 



T(n) 



e-"" a"-^ da (2.2) 







(a-representation). Multiplying two such representations, we have 
1 1 



a1'-'a^^-'daida2. (2.3) 



a^^a^^ T{ni)T{n2) 
Now we proceed to the new variables ai = xa, 0^2 = {1 — x)a, and obtain 

1 

1 T{ni + n2) f a;"i"^(l -a;)"2-irfa; 



afa^^ T{ni)T{n2) J [a^x + 02(1 - x)^ 



hn2 



(2.4) 







This Feynman parametrization is valid not only when Reai_2 > 0, but, by the analytical 
continuation, in all cases when the integral in x is well-defined. 

Combining the denominators in (p.l|) and shifting the integration momentum k 
k — xp, we obtain 

^_r(ni + n2) r x"i-i(l - x)"2-idxd'^fc 

~ r(ni)r(n2) J [-P + x{i - x){-p^) - zO]"^+"^ ■ 
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Now we shall calculate the one-loop massive vacuum diagram, which appear as a sub- 
expression in many calculations. Making the Wick rotation ko = ikEo and going to the 
Euclidean space k"^ = —k\, we have (here 27r'^/^/r((i/2) is d-dimensional full solid angle) 

df^k 271^^/2 r ktUkE , , 

/ I. - (2-6) 



-P + m2 - iO)" r(ci/2) } {k\ + m^)" 

At this point, the meaning of the d-dimensional integral for an arbitrary d becomes com- 
pletely well-defined. Then we use IkEdk^ dk^, and proceed to the dimensionless 
variable z = k\/m^ to obtain 

oo 



T{d/2)' ' J (z + l)^ 



The substitution x = l/{z + l) reduces this integral to the Euler B-function, and we arrive 

at 

. p/t -n^n = .vr^/2 ^(-^f + ") (^2)^/2-» . (2.7) 
^— + 171^ — zO)"- 1 [nj 

We can now resume the calculation of the integral G ( p.5|) . We shall assume that < 0, 
so that production of a pair of on-shell massless particles is not possible. Using { \2.T\ ) with 
— > x(l — a;)(— p2) and n — > ni + n2, we obtain 



1 

^ ^ ^^,/2 r(-^/2 + + 712) ^_^2)d/2-n,-n, fd/2-n,-l _ ^y/2-n,-l 

r(rai)r(n2) y 



This integral is a B-function. Our final result can be written as 



(2.8) 



D, = -{k + p)\ D2 = -k\ (2.9) 
r(-d/2 + ni + n2)r(d/2 - ni)T{d/2 - ria) 



G{ni,n2) 



T{ni)T{n2)T{d - rii - 722) 



where the correct iO terms are assumed in -Di,2- If ^1-1,2 are integer, G{ni, 71,2) is proportional 
to Gi = r(l + e)r2(l — e)/r(l — 2e), the coefficient being a rational function of d. 

Ultraviolet (UV) divergences of the original integral ( p.l| ) are reproduced by the loop- 
momentum integral (|2.?| ), and are given by the poles of r{—d/2 + ni + n2) (with minus sign 
in front of d). Infrared (IR) divergences of the original integral reside in the integral in the 
Feynman parameter x ( p.8|) , and are given by the poles of r{d/2 — ni) and r{d/2 — 712) 
(with plus sign in front of d). 

In many cases, calculation of diagrams in the coordinate space can be simpler than 
in the momentum space. In particular, the one-loop propagator diagram of Fig. ffl in the 
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coordinate space is just the product of two propagators. The massless propagators in the 
p-space and the x-space are related to each other by the Fourier transform: 



d^X = -il-^^-R-^l^ i, ^ 7 TTFT- 2.11 



'-x^ + iO)" r(n) (-p2 - i0)'^/2 

(sanity checks: transform l/(— — zO)" to x-space (|2.1(]| ) and back (p.ll|) , and you get 
l/(— p2 — iO)"; take the complex conjugate of ( 2.1CI|) , rename x p, n (i/2 — n, and 



multiply by {^t^Y i you get ( p.llp ). Multiplying two propagators (|2.10|) with the degrees 
n\ and n^-, we get 

,2(d-n,-n.)^. r(d/2-ni)r(rf/2-n2) 1 

r(ni)r(n2) (_a;2)d-ni-n2 " 

Applying the inverse Fourier transform ( |2.11| ), we reproduce the result ( ^.91 ). 

The one-loop diagram ( |2.9| ) is an analytical function in the complex plane with a 
cut. The cut at p2 > 0, where the real pair production is possible, begins at the branching 
point at the threshold = 0. This cut comes from the factor (— p^)"*^, which is equal to 
(p^)~^ e"""^^ at the upper side of the cut and (p^)"*^ e*'^^ at the lower side, and hence has the 
discontinuity — 2i(p2)-^ sin vre. Therefore, the discontinuity of the one-loop diagram ( |2.9| ) 
at e — is proportional to the residue of G(ni,n2) at e = 0. The discontinuity of the 
diagram with rii = n2 = I can be directly calculated using the Cutkosky rules: draw a cut 
across the loop indicating the real intermediate state, and replace the cut propagators by 
their discontinuities 

^^^-2m5ip'). (2.12) 

3 One-loop HQET propagator diagram 

Now we shall calculate the one-loop HQET propagator diagram with arbitrary degrees of 
the denominators (Fig. 0) 

(-(A; + p)o-^0)"i(-A;2-20)"2 ■ ^^-^^ 

It depends only on a; = po, and not on p. 

We cannot use the ordinary Feynman parametrization ( |2.4| ), because the denominators 
have different dimensionalities. Therefore, we make a different change of variables in the 
double a-parametric integral ( p.3| ): ai = ya, a2 = a, and obtain the HQET Feynman 
parametrization 

oo 

1 T{n,+n2) f y^'-^dy 



a 
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Figure 2: One- loop HQET propagator diagram 



If the denominator ai has dimensionahty of energy, and 02 — of energy squared, then the 
Feynman parameter y has dimensionahty of energy; it runs from to 00. Combining the 
denominators in (|3.1|) and shifting the integration momentum k ^ k — yv/2, we obtain 

r(ni)r(n2) J [-k^ + y{y/A - 07) - zO]"^+"^ ' ^ ' ' 

We must have a definite sign of zO in the combined denominator, therefore, iO terms in 
both denominators must have the same sign. 

We shall assume that < 0, so that production of a pair of on-shell particles is not 
possible. Using (2^) with iv? — > y{y/^ — uj) and n — > rii + ^2, we obtain 

CXD 



r(ni)r(n2) 



/ r /y \-|(i/2-ni-n2 
y"-i[y(|-c.)] (3.4) 



We proceed to the dimensionless variable z = y/{—uj). Then the substitution z/4+1 = 1/x 
reduces this integral to the Euler B-function: 



00 



d/2-n.-l (y_ _ ^ ,j.2n, . .d-n,-2n, ^{'d + + 2^2)1(^/2 - U^) 

^ V4 y " ^ ' r(-rf/2 + ni + n2) 



(3.5) 

Our final result can be written as 



Di = {kv + uj)/uj, D2 = -k^, (3.6) 

r(-c/ + ni + 2^2)1(^/2-^2) 



I{ni,n2) 



r(ni)r(n2) 



If ni,2 are integer, I{ni,n2) is proportional to Ji = r(l + 2e)r(l — e), the coefficient being 
a rational function of d. 

The integral (|2.7|) in d'^k can have UV divergence, and it produces r{—d/2 + rii + 77-2) 
in the numerator of ( |3.4| ). However, it is too optimistic with respect to UV convergence: 
for example, at rii = 2, n2 = 1 it has no pole at d = 4, while the original integral ( p.l|) is 
UV divergent. The HQET Feynman parameter y has dimensionality of energy and runs 
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up to oo. It is natural to expect that this can produce an extra UV divergence. The region 
y ^ oo produces r{—d + ni + 2^2) in the Feynman parametric integral ( p.5|) . The "wrong" 
UV F-function is cancelled with the denominator of (|3.5|) , and rii + 2n2 determines the 
correct UV behaviour. The region |/ — > produces the IR r{d/2 — 71,2). The rule about 
negative/positive sign of d in UV/IR F-functions remains valid. 

The one-loop propagator diagram of Fig. |^ in the coordinate space is just the product 
of two propagators. The HQET propagators in the p-space and the x-space are related to 
each other by the Fourier transform: 

+00 

^ ^ -r-i^(t), (3.7) 



-uo - my 27r F(n) 



-00 
+00 







(-'r'r(» + i)^ (3, 



The first integral is non-zero only at t > 0, when we close the integration contour down- 
wards; for integer n, it follows from the residue at a; = — iO. In the second integral, we 
substitute — > a; + iO for convergence. Multiplying the HQET propagator ( |3.7| ) with the 
degree ni and the massless propagator ( |2.10|) with the degree 77-2, we get 



-2n2^d/2 



rW2-n2).,,.,,+2„2-d 



F(ni)F(n2) 



Applying the inverse Fourier transform (|3.8| ), we reproduce the result ( p.6|) . 

The one-loop diagram ( |3.6| ) is an analytical function in the complex uj plane with a 
cut. The cut at > 0, where the real pair production is possible, begins at the branching 
point at the threshold = 0. This cut comes from the factor (— u;)"^'^, which has the 
discontinuity —2iuj~'^'^ sin 2'7re. Therefore, the discontinuity of the one-loop diagram (|3.6| ) 
at e — s> is proportional to the residue of I{ni,n2) at e = 0. The discontinuity of the 
diagram with ni = n2 = 1 can be directly calculated using the Cutkosky rules ( p.l2|) and 

^ -27fz6{po) . (3.9) 



4 Two-loop massless propagator diagrams 

There is only one generic topology of two-loop massless propagator diagrams. Fig. If 
one of the lines is shrunk into a point, the diagrams of Fig. ^d, c result. If any two adjacent 
lines are shrunk into a point, the diagram contains a no-scale vacuum tadpole, and hence 
vanishes. 
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Figure 3: Two- loop massless propagator diagram 



We write down the diagram of Fig. 13 



n"i n"2 n"3 r)'^4 n^s 

^1 ^2 ^3 ^5 



= -kl , 



Do 



-k^ 



D, = -{kl 



(4.1) 



k,? 



It is symmetric with respect to (1 ^ 2,3 ^ 4), and with respect to (1 ^ 3,2 4). If 
indices of any two adjacent lines are non-positive, the diagram contains a scale-free vacuum 
subdiagram, and hence vanishes. If = 0, this is a product of two one-loop diagrams 
(Fig. |b): 



G{ni, n2, ris, n^, 0) = G(ni, n3)G{n2, n^) . 



(4.2) 



If = (Fig.|c), the (3,5) integral (U) gives G(n3, n5)/(-fc|)"3+"5-rf/2. ^his is combined 
with the denominator 4, and we obtain 



G(0, n2, ris, n^, n^) = G{n3, n^)G{n2, n4 + ^3 + ns - d/2) . 



(4.3) 



Of course, the cases n2 = 0, 71,3 = 0, n4 = follow by the symmetry. 

When all > 0, the problem does not immediately reduce to a repeated use of the 
one-loop formula ( |2.9| ). We shall use a powerful method called integration by parts [Q. 
It is based on the simple observation that any integral of d/dki{- ■ ■) (or d/dk2{---)) 
vanishes (in dimensional regularization no surface terms can appear). From this, we can 
obtain recurrence relations which involve G{ni, n2, n^, n^, n^) with different sets of indices. 
Applying these relations in a carefully chosen order, we can reduce any G(ni, n2, n3, n4, n^) 
to trivial ones, like (|43), (^. 

The differential operator d/dk2 applied to the integrand of (^4.1| ) acts as 



^-.^2{k2 + p) + ^2k2 + ^2{k2-k,] 

ok2 D2 -D4 -D5 



(4.4) 



Applying {d/dk2) ■ k2 to the integrand of ( f4.1| ), we get a vanishing integral. On the other 
hand, from (U), 2^2 ■ ^2 = -2D4, 2(^2 + p) ■ k2 = (V) - ^2 - A, 2{k2 - k,) ■ ^2 = 
D3 — D4 — D^, we see that this differential operator is equivalent to inserting 



d-n2-n,-2n4 + ^(( V) - A) + ^{Ds 

J^2 



Da 
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under the integral sign (here {d/dk2) ■ k2 = d). Taking into account the definition 
we obtain the recurrence relation 

{d-n2-n^- 2n4,)G{ni, n2, na, n^, n^) 

+ n2 [G{ni,n2 + 1, ng, n^, n^) - G{ni, n2 + 1, ng, - 1, ^5)] 

+ ^5 [(7(^1, 77,2, ng - 1, n4, + 1) - n2, ng, n4 - 1, + 1)] = . 

This relation looks lengthy. When using the integration-by-parts method, one has to work 
with a large number of relations of this kind. Therefore, special concise notations were 
invented. Let's introduce the raising and lowering operators 

n2, na, n^, n^) = G{ni ± 1, n2, n^, n^, n^) , (4.5) 

and similar ones for the other indices. Then our recurrence relation can be written in a 
shorter and more easily digestible form 

[d-n2-n5- 2^4 + n22+(l - 4") + n55+(3" - 4")] G = 0. (4.6) 

This is a particular example of the triangle relation. We differentiate in the loop 
momentum running along the triangle 254, and insert the momentum of the line 4 in the 
numerator. The differentiation raises the degree of one of the denominators 2, 5, 4. In the 
case of the line 4, we get — 2D4 in the numerator, giving just — 2n4. In the case of the line 
5, we get the denominator 1^3 of the line attached to the vertex 45 of our triangle, minus 
the denominators D4 and D^. The case of the line 2 is similar; the denominator of the 
line attached to the vertex 24 of our triangle is just — p^, and it does not influence any 
index of G. Of course, there are three more relations obtained from ( [4.6| ) by the symmetry. 
Another useful triangle relation is derived by applying the operator {d/dk2) ■ (^2 — ^i): 

[rf - 712 - 724 - 2n5 + n22+(l" - 5") + n44+(3~ - 5")] G = 0. (4.7) 

One more is obtained by the symmetry. Relations of this kind can be written for any 
diagram having a triangle in it, when at least two vertices of the triangle each have only a 
single line (not belonging to the triangle) attached. 

We can obtain a relation from homogeneity of the integral ( [4.1|) in p. Applying the 
operator p- (d/dp) to the integral (^?T|), we see that it is equivalent to the factor 2{d — J2^i)- 
On the other hand, explicit differentiation of the integrand gives — (ni/L'i)(— + Z^i — 
D3) - (722/1^2) (V + D2- D^). Therefore, 

[2{d -ns-n^- n^) - rii - n2 + rail+(l - 3~) + n22+(l - 4")] / = . (4.8) 



This is nothing but the sum of the {d/dk2) ■ ^2 relation ( |4.6|) and its mirror-symmetric 
[d/dki) ■ ki relation. 



Another interesting relation is obtained by inserting (ki+p)^ into the integrand of (|4.1| ) 
and taking derivative d/dp^ of the integral. On the one hand, the vector integral must be 
proportional to p^, and we can make the substitution 

{ki+p)-p ( Di- p 
ki+p-* 5 p = 1 + 



p2 y _p2 J 2 
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in the integrand. Taking d/dp^ of this vector integral produces ( |4.1| ) with 

(Id- En.) + 

inserted into the integrand. On the other hand, the exphcit differentiation in p gives 
d + ^2(A;i + pf + ^2(^2 + p) ■ ih +P), 2{k2+p)- ih +p)=D,-D,-D2. 
Therefore, we obtain 

[irf + rii - - n4 - + (frf - ^ (1" - 3") + n22+(l- - 5")] G = . (4.9) 

Three more relations follow from the symmetries. 
Expressing (^(ni, ?t,2, ^13, n4, 71,5) from ([4.7D: 

^ 721,^2, ns, n4, rig = ^— ^ ^ '-G , 4.10 

we see that the sum rii + 213 + ng reduces by 1. Therefore, applying ( ^.101 ) sufficiently 
many times, we can reduce an arbitrary G integral with integer indices to a combination 
of integrals with ns = (Fig. ([4. 2D), rii = (Fig. (|4.3D), n3 = (mirror-symmetric 
to the previous case). Of course, if max(n2,n4) < max(ni,n3), it may be more efficient 
to use the relation mirror-symmetric to ( [4.7|) . The relation ( [4.9|) also can be used instead 
of ( |4.7| ). Thus, any integral G(ni, ?t,2, 'n.3, 724, 725) with integer iii can be expressed as a linear 
combination of G\ and (^2, coefficients being rational functions of d. Here the combinations 
of F-functions appearing in n-loop sunset massless diagram are 

F(l - (n l)e) 

Methods of calculation of three-loop massless propagator diagrams are considered in [H. 



5 Two-loop HQET propagator diagrams 

There are two generic topologies of two-loop HQET propagator diagrams. Fig. |^a, b. If 
one of the lines is shrunk into a point, the diagrams of Fig. |^c, d, e result. If any two 
adjacent lines are shrunk into a point, the diagram contains a no-scale vacuum tadpole, 
and hence vanishes. 

We write down the diagram of Fig. 

'^'^^ -^'{-l^f^'-^-^-^^lin,, n2, ns, n4, n,) , 



12 3 4 5 

Di = (ki+p)- v/u , D2 = {k2+p)-v/u, ^ ■ ^ 

D3 = -kl , A = -kl , = -(fci - k2f . 
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5 

b 




Figure 4: Two-loop HQET propagator diagram 

It is symmetric with respect to (1 <-> 3, 2 4). If indices of any two adjacent lines are 
non-positive, the diagram contains a scale-free vacuum subdiagram, and hence vanishes. 
If ns = 0, this is a product of two one-loop diagrams (Fig. |c): 

/(ni, n2, ng, n^, 0) = /(ni, ng) /(ng, Ui) . (5.2) 

If ni = (Fig. |d), the (3, 5) integral (^) gives ^(ns, n5)/(-A;2)"3+"5-d/2. ^^jg jg combined 
with the denominator 4, and we obtain 

J(0, n2, ns, n4, n^) = G{n^, n^)I{n2, n4 + ris ns - d/2) (5.3) 



(and similarly for ^2 = 0). If 71,3 = (Fig. |e), the (1,5) integral ( p.6[ ) gives I{ni,n^) 
/ (— 2ct;)"^+^"^~'^; this is combined with the denominator 2, and we obtain 

/(ni, n2, 0, n4, n^) = I{ni,n^)I{n2 + ni + 2n^ - d, n^) (5.4) 

(and similarly for ^4 = 0). 

When all > 0, we apply integration by parts 0. The differential operator d/dk2 
applied to the integrand of ( p.l| ) acts as 

^ % ^2^2 + ^2(^2 - h) . (5.5) 



dk2 D2 <jJ 

Applying {d/dk2) ■ k2 and {d/dk2) ■ {k2 — ki) to the integrand of ( p.l|) , we get vanishing 
integrals. On the other hand, from ( |5.5|) , k2v/uj = D2 — 1, 2{k2 — ki) ■ k2 = D3 — D4 — D^, 
we see that these differential operators are equivalent to inserting 

1^2 

d-n2-n5-2n4 + — + -j^iDs - D^) , 

n2 ^4 

d - 712 - 714 - 2^5 + — -Dl + -7^{D3 - D5) 
L>2 i^4 

under the integral sign. Therefore, we obtain the triangle relations 

[d-n2-n^- 2^4 + ^22+ + n55+(3" - 4")] / = , (5.6) 
[d-n2-ni- 2n5 + n22+l- + n44+(3- - 5")] / = (5.7) 
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(two more relations are obtained by (1 3, 2 ^ 4)). Similarly, applying the differential 
operator 2uj{d/dk2) ■ v is equivalent to inserting 



-2^ + ^Au\D2 - 1) + ^Auj^D, - D,) , 

and we obtain (taking into account the definition ( p. ID 

[-2n22+ + n44+(2" - 1) + n55+(2" - 1")] / = (5.8) 

(there is also the symmetric relation, of course). 

We can obtain a relation from homogeneity of the integral ( p.l| ) in u. Applying the 
operator u{d/duj) to the integral (|5.1| ) multiplied by (— to')""'!""'^, we see that it is equivalent 
to the factor 2{d — ns — — n^) — rii — n2. On the other hand, explicit differentiation of 
(— u;Di)~"^(— u;i52)~"^ gives —rii/Di — n^jD^- Therefore, 

[2(rf - - n4 - ns) - ni - ^2 + + ?222+] / = . (5.9) 

This is nothing but the sum of the [djdk^) ■ k2 relation ( |5.6D and its mirror-symmetric 
(d/dki) ■ ki relation. 

When trying to find the most useful combinations of recurrence relations, it is conve- 
nient to manipulate these relations in an algebraic way. We shall consider the raising and 
lowering symbols, such operators not commuting with rii, etc. Traditionally, 

these operators are written to the right of rij factors. Of course, any relation may be mul- 
tiplied (from the left) by any Ui factors. We can also shift, say, rii rii ± 1 everywhere 
in the relation. This operation can be represented by multiplication by 1^ from the left 
followed by commuting this operator to the right, if we assume 

l±ni = (ni ± 1)1^. (5.10) 

The most useful combination of recurrence relations for the integral / ( ^.ID is the 



triangle relation (|5.?1 ) minus 1 times the homogeneity relation (|5.9D : 
(i — ni — 712 — n4 — 2n5 + 1 
— (2{d — 77-3 — n4 — n^) — — ^2 + l) 1~ + 77,44^(8^ — 5 

Expressing /(ni, 712, 77,3, 714, 115) from this relation: 



0. 



(5.11) 



/(71i,n2,7l3,7Z4,n5) = ^ — X 

d - rii - n2 - n4 - Zn^ + I (5.12) 
[{2{d -ns-rii- 715) - 7ii - 722 + 1)1" + 7i44+(5~ - 3~)] / , 

we see that the sum rii + 113 + reduces by 1. Therefore, applying ( |5.12D sufficiently 
many times, we can reduce an arbitrary I integral with integer indices to a combination 
of integrals with 7x5 = (Fig. |c, (^), 77i = (Fig. |d, {^), 773 = (Fig. |e, {^). Of 



13 



course, if max (722, ^4) < max(ni,n3), it may be more efficient to use the relation mirror- 
symmetric to Thus, any integral I{ni, ^2, n^, n^, n^) with integer can be expressed 
as a linear combination of and /2, coefficients being rational functions of d. Here the 
combinations of F-functions appearing in ra-loop sunset HQET diagram are 

/„ = r(l + 2ne)r"(l-e). (5.13) 

We write down the diagram of the second topology. Fig. ^b, as 



n'rii r-)7i2 7-)n3 p,n4 j^n^ 
^5 



-7r'^(-2^)2('^-"*-"5) J(ni, n2, ns, n,, n,) 



'5 

Di = {h + p) ■ v/uj , D2 = {k2+ p) ■ v/uj , D3 = {ki + k2 + "[)}■ v/u , 
= -kl , = -kl . 



(5.14) 



It is symmetric with respect to (1 ^ 2, 4 5). If n4 < 0, or < 0, or two adjacent heavy 
indices are non-positive, the diagram vanishes. If = 0, this is a product of two one- loop 
diagrams (Fig. |^c). If ni = or n2 = 0, this is a diagram of Fig. When ^1,2,3 are all 
positive, we can insert 

1 = Di + D2- 

into the integrand of (|5.14]) , and obtain 

J= (1- + 2" - 3-)J. (5.15) 

This reduces ni + n2 + by 1. Therefore, applying ( 5.15| ) sufficiently many times, we 
can reduce an arbitrary J integral with integer indices to a combination of integrals with 
rii = 0, ?T,2 = 0, = 0. 

Methods of calculation of three-loop HQET propagator diagrams are considered in [H. 



6 Renormalization of QCD 

The QCD Lagrangian expressed via the bare (unrenormalized) quantities (denoted by the 
subscript 0) is 

i 

where Dof^qQ = {d^ — igoAQ^t°-)qo, gq is the gauge fixing parameter, c is the ghost field, and 
Dq^iCq = {d^d""^ — (7o/"'"^^pCo. The renormalized quantities are related to the bare ones by 

g^o = Z^'^qi , Aq = Z^l^A , Cq = ^^^c , = Zl/^g , m^o = Z^rrii , = , 

(6.2) 
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where renormalization factors have the minimal structure 



e 4:71 V e / KiirJ 



+ ■■■ (6.3) 



The Lagrangian has dimensionality [L] = d, because the action S = J Ld'^x is exactly 
dimensionless in the space-time with any d. The gluon kinetic term, which is the Qq term 
in jGq^^.Gq^'^ , has the structure (dAo)^; hence, the dimensionality of the gluon field is 
[Aq] = 1 — e. Similarly, from the quark kinetic terms, the dimensionality of the quark fields 
is [qio] = I ~ The covariant derivative Dq^ = — igoAQ^t"" has dimensionality 1, hence 
the dimensionality of the coupling constant is [qq] = e. 

We define to be exactly dimensionless: 

Here n is the renormalization scale, and its rescaling Jt (introduced in the MS renormal- 
ization scheme) allows one to get rid of the Euler constant 7 and log Air in all equations in 
the limit e ^ when they are written via a^. The renormalized QCD coupling satisfies 
the renormalization group equation 

d\ogas , . , . IdlogZa as fas\^ , . 

-r. = -2e-2/3 a, , a, = -— = fj + f]^ I \ +... 6.5 

dlogfj, 2 dlogfi 4:71 \47i/ 

where 

11 4 

Po = yC^ - -^TpUf , (6.6) 

and /3i,2,3 are also known. Here Tp = |, Uf is the number of flavours. All the bare 
quantities, including go, are /i-independent. Substituting the expression for as via and 
Za into (|0|), we find 

2 



l-/5.i|+(^o^-^/3..)(£^) (6.7) 



Hence, Z22 contains no new information: it can be expressed via the one-loop term. 
The bare (unrenormalized) gluon propagator —iDo^yij)) has the structure (Fig. ||) 

+ {-i)Dl^{p)iWP{p){-i)Dl^{p)iW\p){-i)Dl{p) + ■■■ 



(6.8) 



where 



^» = ^(^.^-(l-«o)^) (6.9) 
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is the free gluon propagator, and the gluon self-energy (polarization operator) ill^u{p) is 
the sum of one-particle-irreducible gluon self-energy diagrams (which cannot be separated 
into two disconnected parts by cutting a single gluon line). The series ( |6.8| ) implies the 
equation 

To solve this equation, it is convenient to introduce the inverse tensor A~l of a symmetric 



tensor A^"' satisfying A'^A'^^ = 5"^. If 



. . - ^ Tf 



then 



Using these notations, the equation (|6.10|) can be rewritten as 

DoUp) = (dXUp) - ■ (6.11) 

Due to the Ward identity n^,^(p)p'^ = 0, the gluon self-energy is transverse 

n^.(p) = - p>.PuMp^) ■ (6.12) 

Therefore, the gluon propagator is 



Its longitudinal part gets no corrections. The renormalized gluon propagator is related 
to the bare one by DQ^j_y{p) = ZAD^j_y{p). The gluon field renormalization constant Za is 
constructed to make the transverse part of D^^u finite in the limit e — > 0: 

D,Ap) = Ddp'){9,.-'^)+a^^^. (6.14) 

At the same time, it converts to the renormalized gauge parameter a. This is the reason 
why the gauge parameter is renormalized (|6.2|) by the same constant Za- 



m\/\/\/\y\m -\- •v/w\/\.^^^^\/\/\/\y\» -|- •\/\/\/\/\.^^^^/\/\/\/v^^^^/v\/\/^» -)- 



Figure 5: Structure of diagrams for the gluon propagator 
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In the one-loop approximation (Fig. we can use the methods of Sect. ^ to obtain 



P = TpUf - 



2\-e 



{d-2)Gi 



-4P + 



(47r)'^/2 (d-l){d-3){d-4) 
3d - 2 + (rf - l){2d - 7)(1 - ao) - \{d -l)id- 4)(1 - oq) 

4(rf-2) 



(6.15) 



where Gi is defined in ( [4.11|) . Writing Gq^v{p) = ZAGf^^ip), where the gluon field renor- 
malization constant has the minimal form ( |6.3|) , and requiring that the renormalized gluon 
propagator D^^{p) is finite in the limit e — > 0, we can find Za in the one- loop approxima- 
tion. 



oo 



Figure 6: One- loop gluon self-energy 



It is more convenient to present results for anomalous dimensions instead of renormal- 
ization constants, because anomalous dimensions contain the same information but are 
more compact. An anomalous dimension is defined by 



7(«. 



dlogZ 
dlogfi 



a. 



(6.16) 



If Z is gauge invariant, then, differentiating ( |6.3|) and using ( |6.5|) , we find Zn = — |7o, 
Z22 = §70 (70 + 2/30)5 Z21 = If 7o linearly depends on a (as is the case for Za, for 

example), then Z22 contains an additional term from the /^-dependence of a ( |6.2| ) containing 
the one-loop anomalous dimension of the gluon field: 



, 1 as 



+ 



70(70 + 2f3o) + lAo^a - 27ie 
da 



(6.17) 



Hence, Z22 contains no new information: it can be expressed via the one-loop term. For 
the anomalous dimension of the gluon field itself, we obtain from ( |6.15 ) 



a. 



^A=[{a-f)GA + lT^nf\:g + ... (6.18) 
The bare (unrenormalized) quark propagator iSo{p) has the structure (Fig. |^ 



iSoip) = is\p) + z5°(p)(-«)s(p)z5°(p) + ts%p)i-t)j:ip)is%p){-t)j:{p)is''{p) + ■■■ 

(6.19) 
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where 

S'ip) = = 4±^ (6.20) 

is the free quark propagator, and the quark self-energy (mass operator) — ■iS(p) is the sum 
of one-particle-irreducible quark self-energy diagrams (which cannot be separated into two 
disconnected parts by cutting a single quark hne). The series (|6.19|) implies the equation 

Soip) = S\p) + S\p)i:{p)So{p) (6.21) 

with the solution 

1 1 

= (50)-i(p)-E(p) ^ ^-mo- S(p) • ^^-^^^ 

The quark self-energy has the structure 

j:ip)=pi:v{p')+moEs{p'). (6.23) 
Therefore, the quark propagator is 

So{p) = ^ 2\ 1 — n V f 2^ tl^^ > t 2\\ — = ZgS{p) , (6.24) 

1 - Sy(p^) / - (1 - Sy(]9^)) ^(1 + S5(p^))mo 

and the renormalization constants Zq, Zm are constructed to make Zq(l— Sy) and ZmZq{l+ 
S5) finite. 

•-^ — • + •-^ — W^r^ — • + •-> — ^ — ^ — • + • • • 



Figure 7: Structure of diagrams for the quark propagator 
For a massless quark in the one-loop approximation (Fig. |^), 

Extracting the 1/e pole, we find the anomalous dimension 

75 = laCp-^ + ■ ■ ■ (6.26) 

UV divergences don't depend on masses, therefore, this result is also valid for a massive 
quark. Note that there is no UV divergence in the Landau gauge a = 0. 
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Figure 8: One-loop quark self-energy 



7 Renormalization of HQET 

The HQET Lagrangian expressed via the bare (unrenormalized) quantities is 

L = QqW ■ DoQo + Lqcb , Qo = zU'^Q . (7.1) 



All the renormalization constants in ( |6.2|) are the same as in QCD, where the heavy flavour 
Q is not counted in Uf. In order to flnd the new constant Zg, we need to calculate the 
heavy quark propagator in HQET. 

The bare (unrenormalized) static quark propagator ^5*0(01;) has the structure (Fig. ^) 

iSQ{u}) = iS%uj) + i5°(cu)(-z)S(cj)z5°(w) + iS^{uj){-i)i:{uj)iS%uj){-i)i:{uj)iS%uj) + ■■■ 

(7.2) 

where S^{lj) = l/u is the free HQET propagator, and the static quark self-energy (mass 
operator) —iT,{uj) is the sum of one-particle-irreducible HQET self-energy diagrams (which 
cannot be separated into two disconnected parts by cutting a single heavy quark line). 
Summing this series, we obtain 

So{u;) = 1—. (7.3) 

00 — l^yuj) 



Figure 9: Structure of diagrams for the heavy quark propagator in HQET 
In the one-loop approximation (Fig. |T0|), 

After contraction over the indices, the second term in the brackets contains [kvY = [kv + 
uj — oj)"^. This factor can be replaced by cu^, because all integrals without kv + u m. the 
denominator are scale-free and hence vanish. Using the deflnition (|3.6| ), we get 

= -CF ^^l^^t''^ [2/(1, 1) + i(l - ao)/(l, 2 
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and, finally, 



nco) = c /'^~^±~" A = ao-l- ^ 



(47r)^/2 d-A 



d-3 



(7.4) 



where Ji is defined in ( |5.13|) . 




Figure 10: One-loop heavy quark self-energy in HQET 
Therefore, with the one-loop accuracy, the heavy quark propagator in HQET is 

-SoH-I + Cf^^^I^^^S^. (7.5) 

In the coordinate space. 



So{t) = S\t) 



(7.6) 



Now we re-express 5*0 via the renormalized quantities ag, a (this is trivial at the present 
accuracy), and find Zq of the form ( p. 3D from the requirement that the renormalized 
propagator S{lj) = Zq^Sq{u) is finite in the hmit e — > 0. The result for the anomalous 
dimension is 



7Q = 2(a - 3)Cp^ + 



(7.7) 



Note that there is no UV divergence in the Yennie gauge a = 3. 

Two-loop diagrams for the heavy-quark self-energy are shown in Fig. [Tl|. The first one 
contains the one- loop gluon self-energy ( |6.15| ); it can be easily calculated using (|3.6|), and 
is proportional to I2 ( ^.13|) : 



^?o'(-2^) 



l-4e 



{d-2)h 



{A'kY {d-3){d-Ay{d-Q){2d-7) 



4P. 



(7.^ 



The second one is also recursively one-loop; it is proportional to Cpl2- It is also not difficult 
to understand its dependence on the gauge parameter oq. The exterior loop contains the 
heavy-quark line with a non-integer index n (in fact, n = 1 + 2e). The propagator ( |6.9| ) of 
the exterior gluon produces, as compared to the Feynman gauge oq = 1, the extra factor 



1 + 



1 - ao I{n, 2) - 21{n - 1, 2) + /(n - 2, 2) 
~ KnA) 



1 + i(l - ao)(rf - 3) = -^{d - 3)A . 
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The same is true for the interior gluon hne, with n = 1, as in the one- loop case ( |7.4| ). 
Therefore, the result contains A^: 



Figure 11: Two-loop heavy quark self-energy in HQET 

The colour factor of the third diagram can be easily found using the Cvitanovic algo- 
rithm Pi (Fig. 0). The gluon exchange between two quark lines is replaced by exchange by 
a quark-antiquark pair, minus a colourless exchange which compensates its colour- singlet 
part. Correctness of the coefficients in this identity can be checked by closing the upper 
quark line, and closing it with attaching a gluon (Fig. ^3|). Two example applications of 



this algorithms are shown in Fig. |TJ: the first one is a calculation of Cp, and the second one 

shows that the colour factor of the diagram Fig. |ll|c is Ci;' (~2^) = Cp {Cp — \Ca) ■ The 

colour factor of the three-gluon vertex if"-^'^ is defined as via commutator = if"'^H'^ 

(Fig. |T5|). Therefore, the colour factor of the diagram Fig. |ll]d is the difference of those of 
Fig. 0b and Fig. 0c, \CfCa. 

The diagram Fig. |ll]c, after killing one of the heavy-quark lines (|5.15| ), yields or 12'-, 
it contains for the same reason as Fig. |TT]b: 



S2c(^) = Cf [Cf - \Ca) 



-2uo) 



l-4e 



(47r) 



2/2 



(d-4)2 (ci-4)(2d-7) 



A' 



(7.10) 



The diagram Fig. [TT|d vanishes in the Feynman gauge, because the three-gluon vertex 
vanishes after contraction with three identical vectors v. It also vanishes, if longitudinal 



parts of all three gluon propagators are taken, and hence contains no (1 
result is 



ao) term. The 



S2d(cu) = -CfC 



^7o(-2^) 



l-4e 



(4 



II + 



2(2d-7) 



A(l-ao). 



(7.11) 



Collecting these results together, we obtain the bare heavy-quark propagator with the 
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Figure 12: Cvitanovic algorithm 




Figure 13: Checking the coefficients 




Figure 14: Sample apphcations of Cvitanovic algorithm 



. Y. - . n. - . X. 

-> — < — > — — > — *-> — *^ — — ^ > 

Figure 15: Three- gluon vertex 
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two-loop accuracy 



(47r)'^/2 d-A 



(7.12) 



^?o'(-2^) 



-4e 



\{d-2) 



{d-3){d-Q){2d-7) 



TFnfh + 2A^CFh 



AA 
d-3 



CaII 



+ 



{d-2f{d-h) 
(ci-3)2(ci-6) \{d-?,){2d-7) 



+ {d^ -Ad + 5) A -\{d^-M+ 16) - 3)^2 CaI. 



Now we re-express it via as (|6.4|) and a ( |6.2|) ; in the one- loop term, the as corrections in 
Za (|6I7|), ( |6.6| ) and ( |6.17| ), ( |6.18| ) are necessary. After that, we expand in e, and find 
the renormalization constant Zq of the minimal form (|6.3|) such that ^(ti;) = Zq^Sq{ijj) is 
finite at e — > 0. It must have the form ( |6.17| ). We arrive at the anomalous dimension 



a. 



7Q = 2(a - 3)Cf-^ + Cf 

4:71 



3a2-F24a- 179^ 32^ 
6 3 



Att) 



+ ■ 



(7.13) 



Its difference with the QCD quark field anomalous dimension ( |6.26|) is gauge invariant up 
to two loops: 



a. 



lQ-lq = ~6Cf— - Cf 



127 44 
— Ca- SCf - -^TfUj 



An) 



+ 



(7.14) 



8 Heavy Electron Effective Theory 



Now we make a short digression into the abelian version of HQET — the Heavy Electron 
Effective Theory, an effective field theory of QED describing interaction of a single electron 
with soft photons. It is obtained by setting Cf 1, Ca ^ 0, (^q ^ Cq, a. This 

theory was considered long ago, and is called the Bloch-Nordsieck model. 

Suppose we calculate the one-loop correction to the heavy electron propagator in the 
coordinate space. Let's multiply this correction by itself (Fig. 0). We get an integral in 
ti, t'2 with < ti < ^2 < ^5 < t']^ < t'2 < t. Ordering of primed and non-primed 

integration times can be arbitrary. The integration area is subdivided into 6 regions, 
corresponding to 6 diagrams in Fig. |TB|. This is twice the two-loop correction to the 
propagator. Continuing this drawing exercise, we see that the one-loop correction cubed 
is 3! times the three-loop correction, and so on. Therefore, the exact all-order propagator 
is the exponent of the one-loop term: 



Soit) = S%t)exp 



{it/2) 



2e 



_ {AT^yn 



U) 



In this theory, Za = 1, because there exist no loops which can be inserted into the 
photon propagator. Now we are going to show that = 1, too. To this end, let's 
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t'2 t 














+ 



Figure 16: Exponentiation theorem 

consider the sum of all one-particle-irreducible vertex diagrams, not including the external 
leg propagators — the electron-photon proper vertex. It has the same structure as the 
tree-level term: ieov'^T, F = 1 -|- A, where A is the sum of all unrenormalized diagrams 
starting from one loop. Let's multiply the vertex by the incoming photon momentum q^. 



This product can be simplified by the Ward identity for the electron propagator (Fig. \TT\ ): 



iS%^) ieov ■ q iS\p) = ieo^ [p' ■ v - p ■ v] ^ = leo S\^) - S\p) . (8.2) 

p ■ V p ■ V I J 

In the Figure, the fat dot at the end of the photon line means multiplication by q^; a dot 
near an electron propagator means that its momentum is shifted by q. 

Starting from each diagram for E, we can obtain a set of diagrams for A by inserting 
the external photon vertex into each electron propagator. After multiplying by q^, each 
diagram in this set becomes a difference. All terms cancel each other, except the extreme 
ones (Fig. [1^), and we obtain the Ward identity 

A(...') = -?(4^ „ r(.'..)=g°"'-')-g°"'-). (8.3) 

The vertex function is thus also known to all orders. The charge renormalization constant 
Za is obtained from the requirement that the renormalized vertex function goTZ^jl"^ Zq is 
finite. The factor Zq transforms Sq^ in ( |8.3| ) into 5*^^ and hence makes F finite. Therefore, 
the remaining factor [Z^ZaY^"^ = 1 (this is also true in QED). In the Bloch-Nordsieck 
model, Za = ^ and hence Za = 1. 

Due to the absence of charge and photon field renormalization, we may replace cq — >■ e, 
Oo — > a in the bare propagator ( ^.1|) . It is made finite by the minimal (in the sense of (|6.3|) ) 
renormalization constant, which is just the exponent of the one-loop term 

ZQ = exp[-(a-3)£^] , (8.4) 
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Figure 17: Ward identity for the free electron propagator 
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Figure 18: Ward identity for the electron-photon vertex 
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and the anomalous dimension is exactly equal to the one-loop contribution 

7Q = 2(a-3)£. (8.5) 

Note that the electron propagator is finite to all orders in the Yennie gauge. 

What information useful for the real HQET can be extracted from this simple abelian 
model? First of all, we can obtain the Cp term in (|7.12| ) by the Fourier transformation, 
without explicit calculation. There should be no Cp term in the two-loop field anomalous 
dimension ( [7.13|) . The Ward identity for the gluon propagator is more involved than 
Fig. [1^: in addition to the difference of propagators, it contains ghost terms. The Ward 
identity for the quark-gluon vertex in HQET contains an extra term as compared to (|8.3| ). 
At one loop, it comes from the diagram with three-gluon vertex, and is proportional to 
Ca- Therefore, Z^Za is no longer unity. 



9 1/m corrections to HQET Lagrangian 

As discussed in Sec. |l|, we are interested in processes with characteristic momenta and 
energies oj <^ m. The heavy quark effective theory is constructed to reproduce QCD 5"- 
matrix elements expanded up to some order in uo/m. There is a "folk theorem" that any 
S'-matrix having all the required properties follows from some Lagrangian. Therefore, the 
HQET Lagrangian is constructed as a series in 1/m containing all operators having the 
necessary symmetries, with arbitrary coefficients. These coefficients are tuned to reproduce 
several QCD S'-matrix elements expanded to some degree oi uj/m. We should perform this 
matching for a sufficient number of amplitudes to fix all the coefficients in the Lagrangian. 
After that, we can use this HQET Lagrangian instead of the QCD one for calculating other 
amplitudes. ^ 

The HQET Lagrangian is not unique, because the heavy quark field Q can be redefined. 
Such field transformations can be used to eliminate all time derivatives Dq acting on Q, 
except in the leading term ( |1.2| ) |T0[ . 



The velocity v can be varied by a small amount 6v < uj/m without violating applica- 
bility of HQET or changing its predictions. This reparametrization invariance relates 
terms of different orders in 1/m. 

At the level of 1/m terms, the heavy quark spin symmetry and the superflavour sym- 
metry are violated by interaction of the heavy quark chromomagnetic moment with the 
magnetic component of the gluon field. This leads to hyperfine splittings between states 
which were degenerate in the infinite mass limit (such as B and 5*), as well as to violation 
of leading-order relations among form factors. First, we are going to discuss the simpler 
case of a scalar heavy quark. We shall return to the realistic spin | case at the end of this 
Section. 

The only dimension-5 operator in scalar HQET not containing Dq acting on Q is 
Q^D ^Q. Therefore, the Lagrangian is 

L = Q+iDqQ + ^Q+D^Q + --- (9.1) 
2m 
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The additional term is the heavy quark kinetic energy. This Lagrangian leads to the 

-•2 

dispersion law of a free quark pq = pq — m = Ck^. Therefore, at tree level, Ck = 1. The 
Lagrangian ( |9.1| ) can be rewritten in a covariant form: 

= Qttv ■ DQ, - ^Q+DlQ, + ■■■ (9.2) 
2m 

where D± = D — v{v ■ D). More accurately, this term should be written as j.^ where 
0° = —QI^qD^qQvq is the bare kinetic energy operator. It is related to the renormalized 
one as 0° = Zk(fi)OkilJ^), and hence the term in the Lagrangian is where 

The kinetic energy term gives the new vertices (Fig. ^9|): ^2m^J-' '^2m^o'^"(^' ~'~ P')±^ 
i^gl{m^ + ^^^'')9T (where g^^ = g''" - vf'v''). 

lia jjLa yi) 




p p p p' 

Figure 19: Kinetic-energy vertices 

In Sect. ^ we considered the sum of all one-particle- irreducible heavy-quark self energy 

diagrams —iT,{uo) in the infinite mass limit. Let's denote —ii^Tjk{'^-,P±) the sum of all bare 
one-particle-irreducible self energy diagrams at the order 1/m. Each of these diagrams 
contains a single kinetic-energy vertex (Fig. |T9|). Let's consider the variation of S at 
V ^ V + 5v for an infinitesimal 5v {v ■ 5v = 0). There are two sources of this variation. 
Expansion of heavy-quark propagators l/{p-v + iO) produces insertions ipi ■ 6v into each 
propagator in turn. Variations of quark-gluon vertices produce ig^f^dv^ for each vertex 
in turn. Now let's consider the variation of S^, at p± — ^ p± + 5p^ for an infinitesimal 
6p±. No-gluon kinetic vertices (Fig. produce i—pi ■ 6p±; single-gluon kinetic vertices 

m ■ 



(Fig. |T^) produce i^got"'6p'^; two-gluon kinetic vertices do not change. Therefore, 



This is the Ward identity of reparametrization invariance. Taking into account dT^^/dp'^ = 
2{dEk/dpl)p''^ and dE/dv^" = (rfS/rftu)p^, we obtain 

^ = ^ (9 4) 

dpi duo ■ ^ ^ ^ 

The right-hand side does not depend on p\, and hence 

Sfc(u;,pl) = ^^pi + S,o(o^) . (9.5) 
duj 
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This result can also be understood in a more direct way. The momentum p± flows 
through the heavy quark line. No-gluon kinetic vertices are quadratic in it; one-gluon ver- 
tices are linear; two-gluon vertices are p± independent. The p']_ term comes from diagrams 
with a no-gluon kinetic vertex. Terms linear in p± vanish due to the rotational symmetry. 



The coefficient of p]_ in a no-gluon kinetic vertex is i^. Therefore, the coefficient of pj_ 
in the sum of all diagrams is the sum of the leading-order HQET diagrams with a unit 
operator insertion into each heavy-quark propagator in turn. This sum is just —i^, and 



hence we arrive at ( |9.5|) again. 

As we discussed in the beginning of this Section, coefficients in the HQET Lagrangian 
are obtained by equating on-shell scattering amplitudes in full QCD and in HQET with the 
required accuracy in 1/m. The prerequisite of this matching is the requirement that the 
mass shell itself is the same in both theories, with the accuracy considered. The mass shell 
is defined as position of the pole of the full quark propagator. In QCD it is po = ^/m^ + p'^, 
where m is the on-shell mass (see Sect. |10| for more details). To the first order in 1/m, 
this means ^ = 2m- HQET, the mass shell is zero of the denominator of the bare 
heavy-quark propagator 

Soip) = ^0 r T ■ (9-6) 



UJ 



Sfc^)- 



2m 



-*o dT>(uj) -*9 I \^ 



In Sect. 11 we shall obtain at the two-loop level S(0) = and Sfco(O) = 0. I don't know 
if this is true in higher orders or not. If not, these equalities can be restored by adding 
a residual mass counterterm to the Lagrangian, this does not contradict to any general 
requirements. The mass shell is 



dS(w) 
du 



UJ 



uj=0 



CI 
2m 



dS(cj) 
duj 



p' 



(9.7) 



u=0 



It is correct if C° 



1. The minimal ( |6.3| ) renormalization constant 
has to make Ck{fi) finite; here this means = 1. The kinetic energy operator is not 
renormalized; its anomalous dimension is zero in all orders. The coefficient of the kinetic 
energy operator in the HQET Lagrangian is exactly unity. 



(9.8) 



to all orders in perturbation theory, due to the reparametrization invariance! 

In the case of a spin | heavy quark, there is one more dimension 5 operator, in addition 
to the kinetic energy — chromomagnetic interaction |]T2|, |T3| 



Q+iDoQ + ^Q+D 2g - ^Q+B ■ BQ + 
2m 2m 



or, in covariant notations 



2m Am 



(9.9) 



(9.10) 
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where G^u = gG'^^J:"' and a^" = f[7'^,7'^]- In the v rest frame, only chromomagnetic 

components of G^u contribute, because (T°* sandwiched between Q and Q yields zero. 
Again, it is more accurate to write this term as C^O^ = Cm{lj)Om{lj), where = 
Zm{fJ')Om{fJ') is the bare chromomagnetic operator, and = Z~^{fi)Cm{fi) . The kinetic 
term in (|9.9|) does not violate the heavy quark spin symmetry (because it contains no 
spin matrix between and Q), while the chromomagnetic term violates it, producing 
hyperfine splittings. The chromomagnetic interaction gives the new vertices (Fig. 20): 
^gt^cr'^^qu, ^^g'^\t°' ^'t^](^^^ ■ The arguments leading to ( |9.8| ) remain valid in the spin | 
case. The coefficient of the chromomagnetic interaction Cm(yu) is not related to the lower 
order term in 1/m by the reparametrization invariance, and can only be calculated by 
QCD/HQET matching (Sect. H. 

jjLa jjLa ui) 




Figure 20: Chromomagnetic- interaction vertices 



10 On-shell renormalization of QCD 

Now we are going to discuss calculation of on-shell propagator diagrams of a massive quark 
in QCD. Let's write the one-loop diagram with arbitrary degrees of denominators (Fig. ^1]) 
as 

/ r 2 /^-.^.^ .2 -mn. = ^vr^/^m^-^^"^^"^)M(ni, n,) . (10.1) 

J [m^ — [k + mvY — [—k'^ — zO)"^ 

After the Wick rotation k^ = ikEo, k"^ = —k% and transformation to the dimensionless 
integration momentum K = kE/m, it becomes 

Now let's compare it with the one-loop HQET diagram propagator /(ni, (|3.6|) . In terms 
of the dimensionless Euclidean integration momentum K = kE/{—'2Lj), it has the form 

, = 7r'^/2j(^ ^ ) . (10.3) 

(1 -2iKo)"^(^^)"^ V 1, V ; 

The on-shell integral ( |10.2| ) can be cast into the HQET form ( |10.3| ) using inversion []1^ 
K = K'/K'^. The on-shell denominator K^-2iKo = (1 - 2iK'Q) / K'^ produces the HQET 
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denominator. The integration measure becomes 
d'^K'l {K'^f. The final result is 

M(„,. „.) = /(„,. d - - „.) = r(^- ni -2n.)r(-d/2 + n.+ „,) ^ 

V[ni)V[d -rii- ^2) 

This result can also be obtained using the Feynman parametrization ( |2.4| ). If n\^2 are 
integer, M(ni, 712) is proportional to Mi = r(l + e), the coefficient being a rational function 
of d. 

k 



k + mv 

Figure 21: One-loop on-shell propagator diagram 

The inversion interchanges UV and IR behaviour of an integral. Therefore, the poles 
of r{d — Til — 2^2) are IR divergences (sometimes called on-shell divergences in this case), 
and the poles of V[—d/2 + ni + ^2) are UV divergences. The usual rule about the sign of 
d in r functions applies. 



There are two two generic topologies of two-loop on-shell propagator diagrams. Fig. |22| a 
b. We shall start from the simpler type M, Fig. |2 

d'^ki d'^k2 



-IT 



d^2(d-En,)jVf(^^^^2,n3,n4,n5) 



2 I n ™2 /T I _„,\2 (10.5) 



Di = — {ki + mvY , D2 = — {k2 + mv) 
Ds = -kl , = -kl , = -{ki - k2f . 

Using inversion, we can relate it to the HQET two-loop propagator diagram. The on-shell 
denominators D12 produce the HQET denominators, just as in the one-loop case. The 
denominator becomes (K^ - K2f = iK[ - K'^f / {K'^ K'^) . We obtain [0 

M(ni, n2, ns, n^, n^) = I{ni,n2, d - ni - - n^, d - n2 - - n^, n^) . (10.6) 

However, this relation is not particularly useful for calculation of M(ni, n2, ^3, ^4, ns), 
because this HQET integral contains two non-integer indices. The integrals M can be 
calculated using integration by parts recurrence relations. This is not so easy as in the 
massless case (Sect. ^) and the HQET case (Sect. |). We shall not discuss the algorithm 
here; it can be found in the original literature |T^, |TB[. The conclusion is that any integral 
M(ni, n2, ns, n4, 77,5) with integer indices can be expressed as a linear combination of 
Mf (Fig. p2|c) and M2 (Fig. p2|d), coefficients being rational functions of d. Here the 



combinations of F-functions appearing in n-loop on-shell sunset diagrams with a single 
massive line are 

_ r(l + (n - l)6)r(l + n6)r(l - 2ne)T-{l - e) 
"~ r(l-ne)r(l-(n + l)e) ' ^^''^ 
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c d e 

Figure 22: Two-loop on-shell propagator diagrams 



The type (Fig. |2 



Da 



— {ki + mvY , D2 = rri^ — {k2 + mvY , = rn? — {ki + k2 + mvY , 



(10. 



is more difficult. The integrals N can be expressed |T^, |T^, using integration by parts, 
as linear combinations of (Fig. p2|c), M2 (Fig. p^d), and a single difficult integral 
A/"(l, 1, 1, 0, 0) (Fig. |2^), with rational coefficients. Instead of using A^(l, 1, 1, 0, 0) as a 
basis integral, it is more convenient to use the convergent integral A^(l, 1, 1, 1, 1). It can be 
expressed via 3F2 hypergeometric functions of the unit argument with indices depending 
on e. Several terms of expansion in e are known ]T5| : 

(10.9) 



iV(l,l,l,l,l) = 7r2log2-fC(3) + 0(e) 



Until now, we discussed on-shell propagator diagrams with a single non-zero mass. 
Starting from two loops, there are also diagrams with loops of a different massive quark, 
say, c quark loops in h quark self-energy, or vice versa (Fig. |^). Such diagrams can be 
reduced [T^, using integration by parts, to two trivial integrals (products of one- loop ones) 
and two non-trivial integrals. They are expressed via 3F2 hypergeometric functions of the 
mass ratio squared; their finite terms at e ^ are expressed via dilogarithms of the mass 
ratio. 

The ffist three-loop on-shell calculation, that of the electron anomalous magnetic mo- 
ment in QED, has been completed recently [0. A systematic algorithm for calculation 
of three-loop on-shell propagator diagrams in QED and QCD, using integration by parts. 



was constructed and implemented in . 

The on-shell renormalization scheme is most convenient for calculation of on-shell scat- 



tering amplitudes. The heavy-quark part of the QCD Lagrangian (|6.1|) can be rewritten 



31 



Figure 23: Two-loop on-shell propagator diagram with two masses 



as 

L = QQ{i]po-m)Qo + SmQoQo + --- (10.10) 

where m is the on-shell mass (defined as the position of the pole of the full quark propaga- 
tor), and Sm = m — ttlq is the mass counterterm (mg = Z^m). We shall consider it not as 
a part of the unperturbed Lagrangian, but as a perturbation. It produces the counterterm 
vertex (Fig. ^4]) i 6m. 

— ^ ^(r^ 



Figure 24: Mass counterterm vertex 
The bare heavy-quark self energy can be decomposed as 

E(p) = mSi(p2) + {^-m) S2(p') . (10.11) 
The bare heavy-quark propagator is then 

1 

^''^^^ " (1 - S2(p2)) ^^-rn)+Sm- mEi{p^) ' 
It has the pole at = if 

5m = mEi(m2). (10.13) 

The mass counterterm 6m is determined from this equation, order by order in perturbation 
theory; = 1 — Sm/m = 1 — Si(m^). Then, near the mass shell, Si(p^) — Sm/m = 
S'j^(m^) (p^ — m^) + • ■ ■ , where S'i(p^) = dT^iijP') / dp^ . The bare quark propagator ( |10.12| ) 
becomes 

1 — 2^2\jn^) — Zm^i^\[m^) — 

where dots mean terms which are non-singular at —>■ m^. We define the heavy quark 
field renormalization constant in the on-shell scheme Zq by the requirement that the 
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renormalized quark propagator S{p) (which is related to So{p) by So{p) = ZqS{p)) behaves 

) -2m2S'i(m2)]-\ 



as the free one (|6.20|) near the mass shell. Therefore, Zq = [1 — S2(m^ 

Now let's explicitly calculate and Zq at the one-loop order. It is convenient ||19| 
to introduce the function 



T{t) = — Tr(^ + l)J:{mv{l + t)) = J^iim^) + [^2{m^) + 2m^E[{m^)] t + 0{f 



(10.15) 



then Z^ = l- r(0) and Z^^ = [1 - T'{0)]~\ We have 

d'^k 1 1 



T{t) = -iCFglj 



{2-KYDi{t)D2 4:171 



Tr(^ + l)7^(;^ + ^ + m)7'- 



9iiv + (1 - 



where Di{t) 



m 



{p + k)"^ and D2 = —k"^. While calculating the numerator, we can 



express p ■ k via -Di(t) and D2, and omit terms with -Di(t), because resulting integrals are 
no-scale. Omitting also and higher terms, we obtain 



T(t) = -zCfqI 



d'^k 



2 

D'2 



d-2 
2m2 



(1-t) 



Note that this result is gauge- independent. Now, taking into account Di{t) = Di + {Di — 
D2 — 2rn?')t + (9(t^), we arrive at 



Therefore, 



rf- 3 



:i-t) + oif). 



ryos 



(47r)'^/2 



d-3 



;io.i6) 



The equality Z!^ = Zq is accidental, and does not hold at higher orders. 

On-shell renormalization of QCD at two loops has been performed in [|T5| (see also |]17| 
for exact (i- dimensional contributions of loops of another massive quark). The Olg^) term 
in Sm = —m{Z^ — 1), found in ( |10.16| ), is necessary when calculating OIqq) diagrams 
containing the counterterm vertex (Fig. 24). Three-loop results have been obtained re- 
cently |T^. The on-shell mass is gauge invariant to all orders the quark field renor- 
malization Zq is gauge invariant at two loops |jl5[ but not at three . 



11 On-shell renormalization of HQET 

Now we shall consider on-shell renormalization of the HQET Lagrangian ( |7. 1|) . HQET 
mass shell is = 0. All loop diagrams without massive particles are no-scale and hence 
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Figure 25: Two-loop on-shell HQET propagator diagram 

vanish. Only diagrams with loops of massive quarks (of other flavours) can contribute. 
Such diagrams first appear at two loops (Fig. 

The following method is used for calculation of such diagrams pT] , [2^ . For two vectors 
a and h in d-dimensional Euclidean space, the average over directions of a (or h) 



(a ■ hY = (a262)"/2 



£ cos" e sin'^-^ Ode _T (^) r (f ) 
r,s\n'-'ede - r(i) r(^) 



'11.11 



for even n (positive or negative), and for odd n. In particular, in a c? = 1 dimensional 
space the right-hand side is just (a^fo^)"/^, as expected. We can use this formula for {k ■ f )" 
in Minkowski scalar integrals, because they are calculated via Wick rotation. The HQET 
propagator in Fig. p5] produces just an additional power of A;^, and we are left with the 



vacuum diagram of Fig. This diagram has been calculated in |^ : 



^<i^2{d-ni-n2-n3) 



{w? -kf- 20)"i(m2 -kj- 20)"2 [-{ki - k2y - iOp 
T{-d/2 + ni + n3)T{-d/2 + na + n3)T{d/2 - n3)T{-d + ni + na + ng 
T{ni)T{n2)T{d/2)T{-d + ni + + 2^) 



:il.2) 







Figure 26: Two-loop vacuum diagram 
Using this method, we see that S(0) = 0, and 



dE{uj) 



duj 



-iCpgl 



u)=0 



d'^k v^'v" n^^ik) 
{27rY{k-vy (P)2 



'iCpgo 



-CpTp- 



d<^k 



2(rf- l)(rf-6) 
{d - 2){d - h){d - 7) 



k' 



{27rY \ {k-vy 



n(A;2) 
Jk^ 

(11.3) 
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where Il^i,{k) 



energy, k'^/{k 



{k'^Qtiv — kfj_k^)Il{k'^) is the massive quark contribution to the gluon self- 
)^ = —{d — 2), and the sum is over all massive flavours. We find the 
on-shell HQET quark field renormalization constant Zq = [l — {dT,{u) / du!)^=Q\ from 
the requirement that the renormalized propagator S = Sq/Zq (see ( |9.6| )) behaves as the 
free one near the mass shell: 



-4e 



2{d-l){d-Q) 
{d ~ 2){d ~ 5){d - 7) ' 



;il.4) 



This renormalization constant is not smooth at rrij — > 0. This discontinuity comes from 
infrared gluon momenta, where HQET does not differ from QCD. Therefore, the QCD 
on-shell quark field renormalization constant Zq has the same non-smooth behaviour at 
mi ^ P, p. 

Now we are going to calculate Sfco(O) (see (|9.5| )). At the two-loop level, it is given by 
the diagrams of Fig. ^ (where the second one also the mirror symmetric diagram), with 
zero external residual momentum [u = 0, p± = 0). We obtain 



Sfco(O) = iCpOo 



d'^k n(A;2 



d-2 + 



k' 



(k 



0. 



:ii.5) 




Figure 27: Two-loop diagrams for Sfco(O) 



12 Scattering in external gluonic field in QCD 



Now we are going to perform matching for the scattering amplitudes of an on-shell heavy 
quark in QCD and in HQET, with the linear accuracy in q/m, where q is the momentum 
transfer. 

It is most convenient to calculate scattering amplitudes in an external field using the 
background field method [Q. In the QCD Lagrangian ( |6.1|) , we substitute Aq — > H-Aq, 



where Aq is the external field, and choose the gauge fixing term (^D^Aq) /(2ao), where 
= dfj, — igoAQ^t"-. The ghost term is changed correspondingly: 



L = ^qioiiJPo 



miojqio 



1 
2^ 



{D,A^,Y + {D-4)iD[;4 



'12.11 
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Some vertices with the background field Aq differ from the ordinary ones. In particular, the 
gauge fixing term contains a AqA"^ contribution, altering the three-gluon vertex (Fig. 
to 



9of 



010203 



ki-\ k2 

ao 



^2 



9 



ki 



k2 k3 



9 



M1M2 



appear. 

Hi ai jjL a 



It contains no AqA'^ contribution, so that the four-gluon vertex does not change. The ghost 
term in ( |12.1| ) gives the vertices (Fig. p8|b, c) 

-9^r'''' {ki + k2Y , igir^'-^r-'-^g^-^^- . 

The terms with 1/ao in the three-gluon vertex contain fc^^ or fcg ^ ; when they are multiplied 
by the propagator Do^^u{ki) or Do^^^{k2), correspondingly, they extract the term with ao 
from the propagator (|0|), and no terms with negative powers of the gauge parameter ao 



H2 k2 as " ki k2 ~^ ^ 

a he 

Figure 28: Vertices of interaction with the background field 

The sum of one-particle-irreducible vertex diagrams not including the external propaga- 
tors is the proper vertex igof^T^ij), q), where p is the incoming quark momentum, p' = p + q 
is the outgoing quark momentum, /i and a are the background-field gluon polarization and 
colour indices. The vertex function is r^(p, g) = 7^ + A^{p,q), where A^{p,q) contains 
one-loop (Fig. ^) and higher-loop corrections. 

Now we are going to derive the Ward identity for A^(j9, q)q^. Background-field vertices 




obey simple identities shown in Fig. where the colour structure is singled out as the 
first factor, the fat dot at the end of the background gluon line means multiplication by 
q^, a dot near a propagator means that its momentum is shifted by q (as in Fig. [T7| ). 
Starting from each diagram for S(p), we can obtain a set of diagrams for by attaching 
the background-field gluon to each possible place. For example, starting from the one-loop 
diagram Fig. |^ for S(p), we can attach the external gluon either to the quark line, or to 



the gluon one, obtaining the diagrams Fig. Results of their multiplication by q^ are 



shown in Fig. PT|a, b. The differences in the square brackets are equal. The colour factors 
combine to give the colour factor of Fig. ^ times (Fig. |3T|c), due to the definition of the 
colour factor of the three-gluon vertex (Fig. |15]). Therefore, we have the Ward identity 



A^(p,g)g, = S(p)-S(p + g) or T^{p, q)q, = So\p + q) - So\p) . (12.2 
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It also holds at higher orders of perturbation theory. To verify this, one has to derive 



identities similar to Fig. ^ for other background-field vertices. For an infinitesimal g, we 
have 

A''(p.O) = -f^^^^ or np.O) = ^. (12.3) 



The Ward identities ( |12.2|) , (|12.3| ) are very simple, exactly the same as in QED (or 



in the Heavy Electron Effective Theory, Sect. Multiplying the ordinary three-gluon 
vertex by gives an identity which, in addition to the simple difference of Fig. |30|b, has 
additional ghost terms. Therefore, the Ward identities for the ordinary quark-gluon vertex 
function are more complicated than ( [12. 2| ), ( [12. 3| ). 

The renormalized background field is related to the bare one as Ao = Z]I^A. The renor- 

mahzed matrix element is the proper vertex gof^T^^^p, q) times ^^2'^ : ZqZ ^ Zj gt°-T^^{p, q). 



At g = 0, the factor Zq converts ^ in ( |12.3| ) into S ^, making it finite. Therefore, 



ZoiZa = 1, just as in QED (or in Heavy Electron Effective Theory, Sect. H). In other 
words, QqAq = gA. 

Scattering amplitude of the on-shell quark in an external field is gZQu{p')r^t"'u{p), 
where = p'"^ = rn^, — m)u{p) = 0, (^' — m)u{p') = 0. It is UV finite, but may 
contain IR divergences. It can be expressed via two scalar form factors. For comparing 
with HQET, it is most convenient to use the Dirac and chromomagnetic form factors: 

u{p')TnMp) = u{p') (^F°(g2)i^±|^ + G^(g2)l|^^ ia^(p) ^ (12.4) 

where the renormalized form factors are Fi{q^) = ZgF^i^q'^), Gm{q^) = ZqCl^^q"^). The 
form factors can be singled out by the appropriate projectors. Let's rewrite ( |12.4| ) as 



m(p')W^(p) = uip')J2^i'^t^''^(P) 



Fi = FiV), F, = Gi{q'), rr = ^^i-^, n 



2m Am 
and calculate the traces 

1 
4 

Solving the linear system for Fj, we obtain 
1 



Trr^(|^ + l)7;^(^ + ^/m+l) 



2(^-2) 



xlTrr.(^-,l) l^l±p(. + ^)%^l(^ + ^ + l), (12.5) 



C(g') = - , -Trr.(|^ + i) {^^(v + ^X 

""^^ ' 2(rf-2)4 Mi_^V 2mJ 
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Figure 29: One-loop proper vertex 
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Figure 30: Ward identities for the background- field vertices 
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Figure 31: Ward identity for the background- field vertex function 
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There is no singularity in G^(g^) at 
expanding 



0. We can expand the form factors in 



+ rf ^ + ■ ■ 



sphtting q = [q ■ v)v + q± (where q ■ v = — g^/(2m)) and averaging over directions of q± in 
the {d — l)-dimensional subspace orthogonal to v: 



2m 

Thus we obtain 

F^iq') 4Trr{;(^ + i)tv + o(gV"^^ 



4m^ 



nl3 



- 1-d 



Am? 



Gl{q') 



^ iTrr^(^ + l)(7^-t;, 



:i2.6) 



+ 



(d- l)(d-2)4 

The Dirac form factor at g = is unity, due to the Ward identity ( |12.3|) : 



Fi(0) = Z, 

= Z-[1-T' 



1 + -Tr A^(mt;,0)(|^ + l)v^' 



d 1 



TrS(p)(/; + 1^ 



p=mv 



[12.7) 



(see ( |10.15 )). The total colour charge of the heavy quark is not changed by renormalization. 

The heavy-quark chromomagnetic moment fig = G'm(O) (|12.6| ) can be calculated as 
follows. First, we apply the projector in (|12.6|) to the integrand of each diagram, differ- 
entiating all g-dependent propagators and vertices (Fig. ^), and obtain the bare fi^ via 
scalar integrals (|10.4|) . Then we multiply it by Zq ( |10.16|) . The one-loop result is ||12|| : 



flg=l + 



glm- 



Tie] 



(47r)'^/2 2(c/-3) 



[2{d -A){d- h)CF - {d^ -8d + U)Ca\ + 



12.8) 



Setting go — > cq, Cp — > 1, Ca 0, we reproduce the electron magnetic moment fi 



1 + 



2n 



+ 



in QED. It is convergent. The heavy-quark chromomagnetic moment in 



QCD ( |12.8|) contains an IR divergence with the colour factor Ca- The two-loop correction 
to ( |12.8| ) was calculated in [^, and the effect of another massive flavour (say, c loop 
correction to the b quark chromomagnetic moment) was considered in [jl7|. 



13 Scattering in external gluonic field in HQET 

First, let's consider the leading (zeroth) order in 1/m. Let the sum of one-particle- 
irreducible bare vertex diagrams in the background- field method be igQt"'T^{uj^ q), where u 
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and uj' = uo + q ■ V are the residual energies of the initial and the final quark; T^{uj, q) = 
+ K^iuj^q). There are two vectors in the problem, and q^, and hence the vertex 
function has the structure 



As{uj, uj', q^) = Asiuj', uj, q^) , Ka{uJ, J, c^) = -Aa(cf', uj, q^) . 
It obeys the Ward identity 
l^{uj,q)q, = l,{u:' -u:) + laq^ = ^{lo)-^{lo') or V^'{uj,q)q^ = S^\J) -S^\ 



'13.11 



UJ] 



:i3.2) 



or, for g — > 0, 



(13.3) 



On the mass shell u; = 0, cj' = 0, the renormalized scattering amplitude with the linear 
accuracy in q is 



Z°^(l + A,(0,0,0)K(g)t;'^M,(0) 



ryos 



duj 



u^{q)v^u^{0) = Uy{q)v^Uy{0) 



uj=0 



(13.4) 



The total colour charge of the static quark is not changed by renormalization. 

There are two kinds of 1/m corrections to the scattering amplitude: diagrams with a 
single kinetic vertex, and those with a single chromomagnetic vertex. Let's forget for a 
while that we have obtained the result = 1 ( |9.8| ), and denote the sum of one-particle- 
irreducible bare vertex diagrams in the background-field method containing a single kinetic- 
energy vertex igQt°'-^T'^{p, q), where p = uv + p± and = p + q = u'v + p'^ are the initial 

and final residual momenta of the heavy quark; Vk = {p + p'Yx_ + A^(p, q). Dependence on 
p± and p'^ comes only from the kinetic-energy vertex (Sect. which is at most quadratic 
in them. Therefore, the vertex function has the structure 

A^(p,g) =KU^,u:',q'){p + pr^ 

+ Afci(u;,u;',g^)p^ Afei(u;',u;,g^)p'_^ Afco(c^,c^',g^) (13.5) 
+ Ak3{uj,uj',q'^)p']_-Ak3{^',uj,q'^)p'l + Ak2{uJ,uj',q'^) , 
Aks{uj, uj\ q^) = Aks{uj', q^) , 

Afco(u;, uj\ q^) = Ako{^^', uj, q^) , Ak2{'^, u;', q^) = -Ak2i^', Q^) ■ 

Similarly to Sec. |^, the variation of diagrams for A(p, q) aXv ^ v + 5v is equal to that of 
diagrams for Afc(p, q) at p± ^ p± + 5p^, if 5v = —Sp±. Therefore, the reparametrization 
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mvariance ensures 



dcu 

dAa{uj,uj',q'^) 
did 



[13.6) 



The Ward identity 



ensures, due to fpT 



A^(p,g)g^ = S,(p)-S,(p'; 



(13.7) 



- Aks{uJ, u', q^) + Afci(u;, u\ q^){uj' - uj) + Akz{uJ, J , q^)q^ 
AkQ{uj, J , (^){J -uj)\ A/fc2(t^, J , q^)q^ = T.ko{uj) - T^koi^^') ■ 



duj 



:i3.8) 



The first relation here is, due to (|13.6| ), just the derivative of ( |13.2D in oj. 

On the mass shell a; = 0, tu' = 0, the kinetic-energy correction to the renormalized 
scattering amplitude with the linear accuracy in q is 



^Z-(l + A,(0,0,0)K(g)(p + p')^[«.(0) = ^UM){P + P')>M- 



2m 



;i3.9) 



Let the sum of one-particle-irreducible bare vertex diagrams in the background-field 
method containing a single chromomagnetic vertex be 



igot 



(13.10) 



Tm{uj,uj',q'^) = 1 + Am{u! , w' , q'^) . Reparametrization invariance does not relate Am with 
any vertex function of zeroth order in 1/m. We have no better alternative to a direct 
calculation. In order to obtain the on-shell scattering amplitude with the linear accuracy 
in q, we only need the static quark chromomagnetic moment Jl^ = rm(0,0,0). All loop 

diagrams for A^ vanish, except those with loops of some massive quark. Such diagrams 
first appear at two loops. They can be calculated using the method of Sect. [Tl|. The 
renormalized matrix element Jig = ZqJP^ of the chromomagnetic operator is |22| 



at E 



-4e 



<f -M+ 16 



{d - 2){d - h){d - 1) 



13.111 



where the sum is over all massive fiavours (except the heavy fiavour of our HQET, of 
course) . 
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14 Chromo magnetic interaction 



Now we are ready to compare the on-shell scattering amplitudes in full QCD and in HQET. 
HQET spinors u^ip) are 2-component in the v rest frame: '^Uy{p) = 0. The corresponding 
QCD spinors u{mv +p) are related to them by the Foldy-Wouthuysen transformation p5| 



u(mv +p)= ( 1 + + Oiv^ Im^) \ Uyiv) . 
2m 



(14.1) 



Expanding the QCD scattering amplitude up to linear terms in g/m and re-expressing it 
via HQET spinors, we obtain 



The HQET scattering amplitude with the 1/m accuracy is 



;i4.2) 



ru.(O) 



;i4.3) 



Both scattering amplitudes ( |14.2| ) and (|14.3|) are renormalized and hence UV finite. Both 
may contain IR divergences. By construction, HQET is identical to QCD in the infrared 
region, so that these IR divergences are the same. For example, if there are no other 
massive flavours in the theory, all loop corrections to /x^ vanish because they contain no 
scale. These zero integrals contain UV and IR divergences which cancel. UV divergences 
are removed by and IR ones match those in fig. 

Comparing the coefficients of q^/m, we again see that = 1 and CkifJ^) = 1 (9 
Comparing the coefficients of [4,'y'^]/Tn, we obtain 



f^9 



;i4.4) 



In the one-loop approximation, re-expressing (|12.8|) via as{fi) ( |6.4|) and expanding in e, we 
obtain 



z;;,\fi)CM = i + 



( 

An 



-2Le 



2Cf+\- + 2]Ca 



m 

L = log — 
/i 



;i4.5) 



The minimal (|6.3|) renormalization constant Z^ making Cm finite is 

a. 



Aire 



+ 



and the chromomagnetic interaction constant is 



-CaL + C, + Ca)'^ + 
An 



;i4.6) 



14.7) 
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Therefore, the anomalous dimension of the chromomagnetic operator and Cm{ni) are []12 
^^ = 2Ca^ + ---, C„,{m) = l + 2{Cp + CA)'^^ + --- (14.8 



Two-loop anomalous dimension was calculated in p6[ within HQET, and (a week later) 



in H by QCD/HQET matching: 

7^ = 2Ca^ + ^CA{nCA - UTFUf) (^j + ■ • • (14.9) 
The anomalous dimension vanishes in QED, where Ca = 0. Two-loop correction to Cmifn) 



was found in 22 . It is most convenient to calculate 



4tt 



An 



+ 



because it contains no large logs. Then we can use the renormalization group equation to 
find Cm(/u) at /i <^ m. 

The scattering amplitude (|14.3| ) does not depend on an arbitrary renormalization scale 
/x; the matrix element of the bare chromomagnetic operator is also /i-independent. 
Therefore, 



Zra^{fj)Cm{^j) = COUSt . 

Differentiating this equality in log /i, we obtain the renormalization group equation 



(ilog/i 



7m(as(yU))Cm(yU) 



14.10) 



If L = log m//i is not very large, it is reasonable to find the solution as a series in 
as{m). Re-expressing 



Air 



An 



m 



7m{asW) = lo—r— + 7i 



A-K 



An 



+ 



we obtain the equation 



dL 



+ 



70—; 1 + ^PqJ^—, + 7i 



An 



An 



An 



+ ■ 
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which can be solved order by order in ^^(m): 



47r 



70 I y - /5o ) - (71 + Ci7o) L + C2 



An 



+ ■ ■ 

14.11) 



If j^L ~ 1, we have to solve the renormalization group equation ( |14.10| ) in another 



way. Dividing it by (ilogas(/i)/(ilog/i ( |6.5|) (at e = 0), we obtain 

0. 



d log as 

The solution of this equation is 

CmifJ') = Cm{m) exp 



7n^(as) das 
2f3{as) as 



(14.12) 



We can expand the ratio 7m(as) / P{c(s) in C(s- Integral of the first term gives log ^^(/i) / as{m). 
Integrals of all the other terms are powers of a^, and we may expand the exponent of these 
terms: 



Cm(/i) 



as{m) 



TO 
"2/3o 



„ as{m) /5o7i -/?i7oas(/u) -as(m) 
J- + ^1 — : 777^ : r 



14.13) 



The fractional power of as{fi)/as{m) in ( |14.13 ) contains all leading logs (a^L)" in the 
perturbative series ( |14.11| ); the correction inside the brackets contains the subleading logs. 
We cannot use the C2 term here until we know 72. 

The largest term missing in ( |14.11D is (^)^L^. the largest term missing in ( |14.13| ) is 

C2 (fj)^- Comparing these errors, we can estimate the value of L at which ( |14.13|) becomes 
a better approximation than (|14.11| ). 

The solution ( |14.12| ) of the renormalization group equation depends on m and fi. We 
can rewrite it as a product of a function of m and a function of fi. Let's subtract and add 
^ inside the integral in ( 14.12| ). the difference may be integrated from 0, and we obtain 



70 ^ 70 



(14.14) 



K{jj) = exp 



asi^l) 



7m(tts) 7o \ da. 



2/?(a,) 2(3oJ a 



Poll - |3lloas{^J) 



+ 



CUm) = C^{m)K-\m) = 1 + Ci + 



Poll - Pilo\ Oisirn) 



+ 
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The most obvious effect of the chromomagnetic interaction is the hyperfine B-B* spht- 
ting: 



TTiB* - rriB 



m. 



(14.15) 



where Hq{h) is the matrix element of the chromomagnetic interaction operator. The prod- 
uct Cm{fJ')fJ'G{f^) is, of course, /i-independent, and hence 



70 



(14.16) 



The quantity fiQ is /i-independent, and hence is equal to times some number; we 
obtain 



rriB* - rriB 



mt 



(14.17) 



We can write a similar equation for m/5« — mr). The quantities fiQ in 6-quark and c- 
quark HQET differ by an amount of order (as(mc)/7r)^ due to decoupling of c-quark loops 
(Sect. |15D. Multiplying (|14.17|) by m^. -|- niB = 2mb + 0{l/mb) and dividing by a similar 
D-meson equation, we obtain [EB] 



m 



B* 



m 



B 



m 



D* 



m 



D 



A 



Poll - Pilo\ Ois{mc) - asirub) 



Ait 



+ - 



+ 



MS 



m.b 



In the interval between and rrib, the relevant number of flavours is n 



f 



uIb* — rn?B 

2 2~ 



9/25 



1 - 



7921 asirric) - oisimb) 



3750 



TT 



+ 



TT 



14.18) 



The experimental value of this ratio is 0.89. The leading logarithmic approximation gives 
0.84; the NL correction reduces this result by 9%, giving 0.76. The agreement is quite 
good, taking into account the fact that the l/rric correction may be rather large. 



15 Decoupling of heavy quark loops 

Let's consider QCD with ni light flavours and a single heavy one, say, c. Processes involving 
only light fields with momenta -C nic can be described by an effective field theory — 
QCD with TLi flavours. There are 1/m" suppressed local operators in the Lagrangian, 
which are the remnants of heavy quark loops shrunk to a point. This low-energy effective 
theory is constructed to reproduce S-matrix elements of full QCD expanded to some order 
in pi/mc- Operators of full QCD can be expanded in l/rrtc in operators of the effective 
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theory. Coefficients of these expansions are fixed by matching — equating on-shell matrix 
elements, up to some order in pi/rric. Matching full QCD with the low-energy one is called 
decoupling; it is very clearly presented in |27|, where references to earlier papers can be 
found. 

In particular, the light fields of full QCD are related to those of the low-energy theory 

by 

q. = cM^'Qi A, = a(/i)i/2A;, c = Cc(^)'/v, (15.1) 

up to 1/ rric suppressed terms, where the fields are renormalized at fi in both theories. The 
coupling constant and gauge-fixing parameter in both theories are related by 

^7 = Ca(/i)'/V, a = CA(/i)a' (15.2) 

(we shall see in a moment, why the last coefficient is (a)- It is more convenient to calculate 
the coefficients which relate the bare fields and parameters in both theories. After that, 
it is easy to find the renormalized ones: 

C'W-gS^j- w-'-i^^"- ^"('''-Ir^"- f^^-^' 



IS 



The transverse part of the bare gluon propagator ( |6.13| ) near the mass shell in full QCD 

W) = ^. Zy-j^y (15.4) 

In the low-energy theory, Z'^^ = 1, because all loop diagrams for n'(0) contain no scale. 
Therefore, the matching coefficient in the relation Aq = Ca^o is 

1 

= ^ = ^ _ n(o) • ^1^-^) 

Multiplying by this coefficient also converts Oq into Qq in the longitudinal part of the 
propagator. 

In the full theory, only diagrams with c-quark loops (Fig. ^) contribute. It is convenient 
to extract 11(0) using 



p=0 



In the one-loop approximation (Fig. |32|a), 11(0) is a combination of one-loop vacuum inte- 
grals (^.71); we obtain 
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Therefore, 



The ratio of the renormahzation constants ( |6.17| ) for the gluon field is 

^ 1. as A as 

where A7^o = is the contribution of a single flavour to the one-loop anomalous 
dimension ( |6.18|) of the gluon field. Finally, we obtain the renormalized decoupling constant 

aW = i + |r,^iog!^i + .... (15.7) 

It is most natural to calculate Ci{mc) which contain no large logs. The dependence on /i 
can be easily restored using the renormahzation group. It is not too difficult to calculate 
the two-loop diagrams Fig. (|32Db, c. They reduce to combinations of two-loop massive 
vacuum integrals (|11.2|) , and are proportional to r^(e). Using also the two- loop anomalous 
dimension of the gluon field, we can obtain 



GK) = 1--TH4C^-Ca)(^^) +■■■ (15.8) 






a b c 

Figure 32: Massive loops in gluon self-energy 

Similarly, for the light quark field (see (|6.24|) ) 

1 



C 

So 



'70s 

y/os 



l-Sy(O) 



;i5.9) 



The only two-loop diagram contributing is Fig. |33|. Expanding the quark self-energy ( |6.23| ) 
at rriQ = in p up to linear terms, we reduce it to the integrals (|11.2|) , similarly to ( |11.3| ) 



Sy(0) = -tCpg'o 



d 



{27lY (P)2 



n r . 1 g^rn-^' , 2(rf - l)(d - 4)(d - 6) 



(15.10) 



(4vr) 



d{d - 2){d - b){d - 1) ■ 
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The ratio of the renormahzation constants ( |6.17| ) for the quark field is 



§ = 1 + ^ (7,oA/?o + \m 



AO—; — 
da 



Aire) 



+ 



a. 



where A/3o = ^7^0 = and Ajgi = —ACpTp are the single-flavour contributions. 

Finally, we obtain |^ 



2 



C,K) = 1 - ^C^T^ ) +... (15.11) 




Figure 33: Two-loop on-shell massless quark self energy 

The decoupling relation for the coupling constant can be derived from considering any 
vertex in the theory. For example, in |^ it was obtained from the ghost-ghost-gluon vertex. 



Let g^Vc be the sum of bare one-particle-irreducible ghost-ghost-gluon vertex diagrams, not 
including the external propagators. For this vertex, we have 



because after multiplication by the three propagators it becomes the Green function of cq, 
Co, Aq. It is most convenient to nullify all the external momenta, after a single differentia- 
tion in the outgoing ghost momentum. In the low-energy theory, F'^ = 1, if we single out 
the colour factor. In full QCD, only diagrams with c-quark loops contribute. They first 
appear at two loops, and there are three of them. We have 

1 7' 

Calculating C,c and Fc, we can obtain the famous result (see [p7[] ) 

C.W^l-(ic.-Hc,)(^^)%... (15.12) 

It means that at ;U > mc, we have a<i(/i) whose running is given by the j3{as) function with 
rii + l flavours; at /i < mc, we have a^(/i) whose running is given by (3'{a'g) with rii flavours; 
at /i = rrici 
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If we consider 6-quark HQET instead of QCD, nothing changes. When all charac- 
teristic (residual) momenta become much less than m^, c-quark loops shrink to a point. 
From ( |11.4| ), we obtain 



r 



Z'^^ ' ' (47r)'^ ^ '{d-2){d-h){d-7) 

The ratio of renormalization constants is, from ( |7.13| ), 



(15.13) 



Finally, we obtain |28 



CQ(m,) = l + -C^T^( ^— 



+ 



(15.14) 



Finally, we shall discuss decoupling of c-quark loops in the 6-quark chromomagnetic 
interaction constant. Both HQET with c-loops and the effective low-energy theory must 
give identical on-shell scattering amplitudes, up to corrections suppressed by powers of 
l/m^. Therefore, from ( |14.4 ) we obtain 



(15.15) 



Here /i' = 1, and /i is given by ( p.3.11| ), with just the c-quark contribution. The ratio of 
the renormalization constants Z„i (|6.171) is 



m. ^ / 



A-Kt) 



(15.16) 



where A/3o = —\Tp and A71 = — | ■ ISCyiTp are the single-flavour contributions. Re- 
expressing ( |13.11D via as(mc) and expanding in e, we see that singular terms cancel, and 



27 V 47r 



+ 



(15.17) 



Therefore, when crossing /i = nic, we have to adjust Cm{lj) according to (|15.17|) , and at 
II < TTic the renormalization-group running is driven by the /3 function and the anomalous 
dimension with nj = 3. 



49 



16 Conclusion 



Part 2 of these lectures will deal with heavy-light quark currents, and with renormalons 
in HQET. More physics applications of HQET can be found in other reviews. Here is a 
partial hst, in no particular order: |^-^ . 



I am grateful to K. G. Chetyrkin, who proposed the idea to give these lectures; to 
Th. Mannel, J. H. Kiihn, T. van Ritbergen for useful discussions during my stay in Karl- 
sruhe; to D. J. Broadhurst, M. Neubert, A. Czarnecki, A. I. Davydychev, O. I. Yakovlev for 
collaboration on various HQET projects; to D. V. Shirkov, D. I. Kazakov, S. V. Mikhailov, 
A. P. Bakulev for organizing the very interesting Calc-2000 school/ workshop in Dubna. 
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